System size dependence of strange particle yields 
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Abstract. Yields and spectra of strange hadrons (K+, K~, (j), A and A) as well as 
of charged pions were measured in near central C+C and Si+Si collisions at 158 A 
GeV beam energy with the NA49 detector. Together with earlier data for p+p, S+S 
and Pb+Pb reactions the system size dependence can be studied. Relative strangeness 
production rises fast and saturates at about 60 participating nuclcons; the net hyperon 
spectra show an increasing shift towards midrapidity for larger colliding nuclei. An 
interpretation based on the formation of coherent systems of increasing volume is 
proposed. The transverse mass spectra can be described by a blast wave ansatz. 
Increasing flow velocity is accompanied by decreasing temperatures for both kinetic 
and chemical freeze out. The increasing gap between inelastic and elastic decoupling 
leaves space for rescattering. 

PACS numbers: 25.75.-q 



1. Introduction 

The production of strange particles in heavy ion collisions is studied in many experiments 
since it had been proposed as a signature for the transition to a deconfined state of 
strongly interacting matter [1]. Although enhanced strangeness production in A+A 
collisions relative to elementary reactions is observed over a wide range of cm. energy 
[2], the question about its origin is still unsolved. This motivated the systematic study 
of symmetric collisions of nuclei with increasing mass number A. 

The presented data were recorded with the NA49 hadron spectrometer [3] at the 
CERN SPS. The analysis of meson (7r='=, K='=, 0) and A hyperon production is described 
in [4], the preliminary results for A baryons were obtained analogously. 

2. Results 

The yields of all particles under study were measured in {y,pT) bins, the transverse mass 
(rriT — ml + p^) spectra at midrapidity are shown in figure 1 (left). Fits with the 
thermal ansatz 
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Figure 1. Left: transverse mass spectra of various particles, measured at midrapidity 
in C+C (■) and Si+Si (• ) reactions. Shown are the exponential fits with equation 
1. Right: rapidity spectra of A (top) and A hyperons (bottom, preliminary). Only 
statistical errors are included. 



were used to extrapolate to px regions not covered by the measurement. The rapidity 
densities dn/dy of the hyperons are shown in figure 1 (right), while the mesons 
are presented in [4]. For the extrapolation to the very forward region the pion 
rapidity distributions were approximated by a superposition of two Gaussians displaced 
symmetrically around midrapidity, whereas single Gaussians were used for the other 
mesons and A hyperons. To obtain the yields of A hyperons their rapidity distributions 
were extrapolated by shapes adopted from p+p and S+S measurements. The mean of 
both approximations was taken for C+C, the shape from S+S was used for Si+Si (lines 
in figure 1). 

3. Discussion 

3.1. Strange particle production 

The strange particle yield per pion is shown in figure 2 together with data from p+p 
and Pb+Pb collisions [5] measured by the NA49 Collaboration and S+S results [6] 
from NA35. All presented data on pions and hyperons are corrected for feed down 
contaminations. The mean number of nucleons in the geometrical overlap area (Nw) 
serves as a measure of the system size. The ratio of strange hadrons to pions exhibits 
for all considered particles a fast rise in small systems and reaches the level of Pb+Pb 
interactions at about 60 wounded nucleons. 

Empirical scaling parameters as e.g. the nucleon or collision density [7] suggest that 
the density reached in A+A interactions has an impact on the strangeness production. 
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Figure 2. Left: multiplicity of hyperons, charged kaons and mesons per pion, 
(tt) = ^{{'^^) + ))• Right: A density at midrapidity per wounded nucleon in p+p 
(▼), p+Be (0), C+C (■), Si+Si (• ), minimum bias (a) and central Pb+Pb (a) 
collisions at 158 A GeV as well as in S+S reactions (if) at 200 A GeV . The error bars 
represent the squared sum of statistical and systematic errors. 



At sufficiently high density one might expect that the subsequent colhsions do not occur 
independently anymore, but that connected domains of overlapping resonances or strings 
are created. These volumes can become quantum-mechanically coherent and decay as 
objects which can be considerably larger than those created in isolated p+p interactions. 

The volume dependence of strangeness production is described in statistical models 
by the transition from canonical to grand canonical ensembles. The calculation by 
Tounsi and Redlich [8] agrees quahtatively with the measurements presented in figure 2, 
but quantitatively the saturation level is reached for a distinctly smaller number of 
participants in the model. Theory and experiment might be reconciled by two plausible 
modifications: assuming that only parts of the reaction volume that is spread over about 
3 units of rapidity [9] are coherent or that the linear relationship between volume and 
number of participants used to link model and measurement has to be refined. 

3.2. Chemical freeze out conditions 

The measured particle yields can be reproduced by the statistical hadronisation model 
of Becattini [10], see figure 3. The strangeness production in the corresponding 
equilibrated resonance gas is suppressed (75 < 1). This suppression diminishes with 
increasing system size (figure 4), following the observed strangeness enhancement. The 
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Figure 3. Comparison of the measured preliminary particle yields in C+C (left) and 
Si+Si collisions (right) to calculations with the statistical hadronisation model [10]. 



baryochemical potential fiB does not depend significantly on the system size, while the 
decoupling temperature T^h is higher in smaller reaction volumes. This would leave less 
room for inelastic rescattering in the small systems, if the varying Tch is not entirely 
provoked by the changing 75 as suggested in [14]. The freeze out points of C+C and 
Si+Si reactions appear in the QCD phase diagram in the vicinity of the phase boundary 
as seen in figure 4. 

3.3. Absorption of A hyperons 

The increasing absorption of anti-baryons (p and d) with increasing centrality of Pb+Pb 
reactions at 158 A GeV was observed in a forward rapidity window (y = 3.7) by the 
NA52 Collaboration [23] and explained with an increasing baryochemical potential. In 
contrast to that, the (A) yield per pion in central Pb+Pb collisions is not significantly 
smaller than in Si+Si reactions (figure 2 left), while the midrapidity densities are on the 
same level. Moreover the measurement by the NA57 Collaboration [24] indicates even 
a further increase from central Si+Si to Pb+Pb data at different centralities (figure 2 
right). This is in accord with the almost constant fis shown in figure 4. 

3.4. Particle ratios 

Both, the (K~)/(K+) as well as the (A) /(A) ratio (figure 5) of the total yields exhibit 
no significant dependence on the system size. Due to that, the strange hadron potential 
/is calculated (with the strange quark potential fig) as [25] 
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Figure 4. Left: chemical freeze out conditions for C+C, Si+Si, Pb+Pb at 158 A 
GeV from fits with the statistical hadronisation model [10, 11]. The curves are shown 
to guide the eye and represent a functional form f{{Npart)) = a + b ■ exp (c • (Npart))- 
Middle: freeze out conditions for p+p, C+C, Si+Si, Pb+Pb at 158 A GeV. The 
temperatures at chemical decoupling are from fits with the statistical hadronisation 
model [10, 11]. The kinetic freeze out temperature and the transverse flow velocity 
Pt are from blast wave fits to the transverse mass spectra, equation 3. The values for 
p+p and Pb+Pb interaction were taken from [12]. The critical temperature Tc is from 
lattice QCD calculations [13]. 

Right: QCD phase diagram. The phase boundary between deconfined and hadronic 
matter was estimated with lattice calculations by the Bielefeld [15] as well as by the 
Budapest group [16]. The parametrisation of the chemical freeze out points as a curve 
of constant energy per hadron of (E) / (N) = 1 GeV is taken from [17]. The C+C and 
Si+Si data [10] lie between the Si+Au measurements at the AGS at 14.6 A GeV and 
the S+S data at 200 A GeV. Model calculations with 7s = 1 are from [18], those with 
free 75 are from [19]. The freeze out points determined from the fits to the Au+Au and 
Pb+Pb data with 75 = 1 from SIS [20], AGS [18] and SPS [21] up to RHIC energies 
[22] agree with those from the fits with free 75 for the AGS and SPS [11]. 




Figure 5. (A) / (A) (left) and (K^) / (K+) ratio (middle) of the total yields measured in 
p+p (▼), C+C (■), Si+Si (• ), minimum bias (a) and central Pb+Pb (a) collisions 
at 158 A GeV as well as in S+S reactions {ir) at 200 A GeV. The error bars represent 
the squared sum of statistical and systematic errors. The strange hadron potential 113 
(right) is calculated from these ratios. 
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Figure 6. Left: transverse mass spectra of various particles from C+C (left part) and 
Si+Si collisions (right part) measured at midrapidity. Also shown are the fits with the 
blast wave model (equation 3). 

Middle: Inverse slope parameters of A (■) and A (• ) hypcrons (equation 1) in p+p, 
C+C, Si+Si, S+S and Pb+Pb collisions. 

Right: net hyperon rapidity density normalised to the number of wounded nucleons 
{Nw) measured in p+p (t), C+C (■), Si+Si (• ) and Pb+Pb (A) collisions at 158 A 
GeV as well as in S+S reactions (if) at 200 A GeV. Only statistical errors are included. 

is constant as well. In spite of zero net strangeness the potential is not vanishing, it 
amounts to about 60 MeV (figure 5). This is in agreement with the statistical model 
prediction [26] for the fitted values of T^h and /i^. 

3.5. Net hyperon density 

The A/A ratio at midrapidity (not shown) is (in contrast to the fiat (A) /(A) ratio) 
steeply decreasing with increasing system size for small reaction volumes, followed by a 
saturation above about 60 wounded nucleons. The ratio refiects directly the changing 
A rapidity distribution while the shape of the A hyperon spectra remains the same. The 
difference of the two distributions, the net hyperon density, is shown in figure 6. The 
fiattening of the (A — A) rapidity spectra with increasing system size can be understood 
in terms of stronger stopping due to an increasing number of collisions. This leads to 
a successive shift of the incoming baryons from beam towards mid rapidity. Thereby 
the energy per nucleon deposited in the fireball is increasing and more energy for the 
excitation of resonances or strings is provided. 

3.6. Transverse mass spectra 

The inverse slope parameters T from exponential fits to the my spectra with equation 1 
are systematically higher for A than for A hyperons as can be seen from figure 6. In small 
systems this may be traced back to the anti-correlation of the (higher A) production 
threshold on one hand and the kinetic energy of the created particles on the other hand 
due to energy conservation [27]. The cause in the large systems is not clear. 
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The slopes of both, A and A hyperons, fit into the picture of growing radial flow 
with increasing system size. Moreover the uit spectra can be approximated with a 
hydrodynamical model; here a simplified version from Schnedermann et al. [28] with 
constant flow velocity /?t is utilised to fit all particles except the pions simultaneously 
with the function 

dn /pr-sinhp\ /mr-coshpx 

oc mT-hi ■ Ki 7f , (3) 



' • "^1 U \ rri I 1 \ rri 

rriT ■ drriT \ I / \ I 

where Iq and Ki arc modified Bessel functions and p includes Pt- P = atanh Pt- The 
fit to the C+C and Si+Si data is shown in figure 6, while the system size dependence 
of the kinetic freeze out temperature and flow velocity Pt is displayed in figure 4. 
The p+p and Pb+Pb spectra were fitted with the same formula. In p+p reactions 
Pt is compatible with zero; in larger systems Pt steeply increases. The increasing 
radial flow is accompanied by a decreasing kinetic freeze out temperature; T^in drops 
by 30 MeV between C+C and Pb+Pb collisions. Thereby the gap between chemical 
and thermal decoupling temperatures A T increases with size of the colliding nuclei, so 
more rescattering is expected in the larger systems. 

The largest A T is observed in p+p interactions, where neither rescattering nor flow 
is expected. This questions the mechanism described above as the only source of the 
effect. The broadening of pr spectra is already seen in p+A data and multiple scattering 
is given there as an alternative explanation [29]. Furthermore there are models which 
describe the measured yields and spectra reasonably well by a single freeze out, e.g. 
[30]. 



4. Summary 

Yields and spectra of strange hadrons and charged pions in p+p, C+C, Si+Si, S+S and 
Pb+Pb coUisions were presented. The fast rise of the relative strangeness production, 
followed by a saturation above about 60 participating nucleons, together with the 
increasing shift of net hyperons towards mid rapidity with increasing size of the colliding 
nuclei suggests the creation of coherent domains as a possible interpretation. The 
hadronisation of the small systems occurs in the vicinity of the phase boundary. No 
significant absorption of A hyperons is seen. The transverse mass spectra can be 
described by a blast wave ansatz. Fits to the data indicate increasing flow velocity 
accompanied by decreasing temperatures for both kinetic and chemical freeze out. The 
increasing gap between inelastic and elastic decoupling in A+A collisions leaves space 
for rescattering. 
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